AHP中互反判断矩阵一致性调整的新简易方法及交互式实现

提出Satty′s互反判断矩阵一致性调整的一种新的简易方法,并借助EXCEL电子表格的动态引用功能,实现了判断者在建立判断矩阵时可即时地进行一致性的互动式调整;不但改善了现有不合理的强制性调整方法,而且提供了一种有效辅助判断矩阵建立的途径.其中,借鉴欧式距离原理构造了"行和列偏离程度系数"来衡量一致性优化调整过程有效性,同时实现了权重向量相对于互反判断矩阵调整的灵敏度分析.最后通过一个具体实例演示方法的操作过程,并检验了其有效性.     

1,3,3-三硝基氮杂环丁烷的热安全性

借助不同加热速率(β)的非等温DSC曲线离开基线的初始温度(T₀)、onset温度(T_e)和峰顶温度(T_p), Kissinger法和Ozawa法求得的热分解反应的表观活化能(E_k和E_O)和指前因子(A_k), Hu-Zhao-Gao方程ln β_i＝ln[A₀/(b_{e0 or p0}G(α))]＋   b_{e0 or p0}T_{ei or pi}求得的b_{e0 or p0}, Zhao-Hu-Gao方程ln β_i＝ln[A₀/((a_{e0 or p0}＋1)G(α))]＋(a_{e0 or p0}＋1) ln T_{ei or pi}求得的a_{e0 or p0}, 微热量法确定的比热容(C_p), 以及密度(ρ)、热导率(λ)和分解热(Q_d, 取爆热之半)数据, Zhang-Hu-Xie-Li公式、Hu-Yang-Liang-Xie公式、Hu-Zhao-Gao公式、Zhao-Hu-Gao公式、Smith方程、Friedman公式和Bruckman-Guillet公式, 计算了TNAZ在β→0时的T₀, T_e和T_p值(T₀₀, T_eo和T_p0)、热爆炸临界温度(T_be和T_bp)、绝热至爆时间(t_Tlad)、撞击感度50%落高(H₅₀)和热点起爆临界温度(T_cr), 得到了评价TNAZ热安全性的结果: T_SADT＝T_e0＝485.81 K, T_p0＝497.38 K, T_beo＝499.50 K, T_bp0＝513.45 K, t_Tlad＝8.90 s (n＝0), t_Tlad＝8.96 s (n＝1), t_Tlad＝9.01 s (n＝2), H₅₀＝28.88 cm, T_cr＝641.46 K (T_room＝293.15 K), T_cr＝658.89 K (T_room＝300 K), 表明: (1) TNAZ对热是稳定的; (2)撞击感度好于环三亚甲基三硝胺(RDX); (3)热点起爆临界温度高于RDX, 而界于1,3,5-三氨基-2,4,6-三硝基苯(TATB)和六硝基茋(HNS)之间.     

非线性最小二乘跟踪的对偶变量变分方法

最小二乘跟踪方法是近几年提出的一种计算动力系统跟踪轨迹的方法.基于最小二乘跟踪的灵敏度分析算法可以有效避免传统的非线性系统灵敏度分析方法中的病态初值问题,因此其在混沌系统灵敏度分析方面有着重要的应用.针对非线性的最小二乘跟踪问题,首先将其重新描述为带有约束的非线性最优控制问题,引入协态变量并将系统的哈密顿函数表示为关于状态变量和协态变量的函数.然后将目标函数的积分时间离散化,根据对偶变量变分原理,以离散区间两端的状态变量作为独立变量,用Lagrange插值多项式近似离散区间内的状态变量和协态变量,进而将非线性最优控制问题转化为求解非线性方程组问题.这种算法无需对原问题做线性化处理,避免了复杂的线性化过程以及可能因此造成的误差,同时为求解非线性最小二乘跟踪问题提供了新的思路.根据最小二乘方法可以得到两条设计参数有微小变化的状态轨迹,基于这两条状态轨迹可进一步计算出系统关于设计参数的灵敏度,范德波振子作为数值算例验证了该方法在求解最小二乘跟踪问题以及计算非线性系统灵敏度时的有效性.      

高超声速舵面热气动弹性不确定性及全局灵敏度分析

鉴于高超声速中气动热预测的不确定性影响热气动弹性分析的可靠性,提出一种温度分布参数化模型,基于此模型,对高超声速舵面热气动弹性中气动热的不确定性及全局灵敏度进行分析,分析方法:求解N-S方程得到物面的温度分布,对此温度分布进行参数化,分别采用蒙特卡罗模拟(Monte Carlo simulation,MCS)方法和稀疏网格数值积分(spare grid numerical integration,SGNI)方法生成不确定性及全局灵敏度分析所需样本,对所有样本都进行热气动弹性分析,热气动弹性分析过程为:由样本得到温度分布,基于此温度分布,考虑热应力和材料属性的影响,对结构进行模态分析,将结构模态插值到气动网格,采用基于CFD的当地流活塞理论进行了气动弹性分析.分别在两种飞行状态下进行分析,计算结果表明:(1) M=5,H=15 km,结构固有频率和颤振分析结果的变异系数约为5.83%;(2) M=6,H=15 km,结构和颤振分析结果的变异系数约为8.84%.两种状态下,两个不确定参数的全局灵敏度都在50%左右,两者耦合作用很小,约为0.与MCS方法相比,SGNI方法显著的提高了不确定性分析效率.      

Two forms of continuum shape sensitivity method for fluid–structure interaction problems

Continuum Sensitivity Equation (CSE) methods for deriving and computing derivatives with respect to shape design variables are developed in two forms and compared in their application to fluid–structure interaction (FSI) problems. The local derivative form poses the CSEs in terms of the partial derivatives of the state variables with respect to shape parameters, while the CSEs in total derivative form are posed in terms of the total derivative, also known as the material or substantial derivative. In the literature CSEs are often posed in local form for fluids and total form for solids. The two forms are compared here for the purpose of applying a single form to both fluid and structure domains. The local form, also known as the boundary velocity method, requires design velocity only at the boundaries and interfaces of the domains to pose the CSEs. In contrast, the total form, also known as the domain velocity method, requires the design velocity in the whole domain. The local form requires higher-order spatial derivatives of the analysis solution than the total form, which affects the accuracy of its results. Higher order p-elements are shown to be a remedy to the inaccuracy of local form CSE seen in the literature for finite element solutions. The practicality, accuracy, and efficiency of these two CSE forms are compared based on the implementation and computed derivatives for three examples: a linear Timoshenko beam subject to a tip force, fluid flow around an airfoil, and an airfoil attached to a nonlinear joined beam subject to a gust load.     

基于荷载时程分析法的钢板混凝土结构墙的抗冲击性能敏感性分析

选取由表面钢板、拉结筋、剪力钉及混凝土组成的钢板混凝土结构墙为研究对象,运用经典显式非线性动力分析软件ANSYS/LS-DYNA,基于荷载时程分析法进行了一系列影响钢板混凝土结构墙抗冲击性能的参数敏感性分析。这些参数包括:墙体厚度、钢板厚度、拉结筋直径与间距等。分析结果表明,以上参数均会影响墙体的抗冲击性能,尤其是墙体与钢板的厚度以及拉结筋的间距。本文的研究工作对于核电厂核岛厂房钢板混凝土结构外墙的设计具有一定的指导与参考意义。     

Sensitivity analysis of pull-in voltage for RF MEMS switch based on modified couple stress theory

An approximate analytical model for calculating the pull-in voltage of a stepped cantilever-type radio frequency(RF) micro electro-mechanical system(MEMS) switch is developed based on the Euler-Bernoulli beam and a modified couple stress theory, and is validated by comparison with the finite element results. The sensitivity functions of the pull-in voltage to the designed parameters are derived based on the proposed model. The sensitivity investigation shows that the pull-in voltage sensitivities increase/decrease nonlinearly with the increases in the designed parameters. For the stepped cantilever beam, there exists a nonzero optimal dimensionless length ratio, where the pull-in voltage is insensitive. The optimal value of the dimensionless length ratio only depends on the dimensionless width ratio, and can be obtained by solving a nonlinear equation. The determination of the designed parameters is discussed, and some recommendations are made for the RF MEMS switch optimization.     

爆轰流体力学模型敏感度分析与模型确认

验证、确认与不确定度量化(VVUQ)是评估物理模型可信度和量化复杂工程数值模拟结果置信度的系统方法.验证是要回答数值模拟程序是否正确求解了物理模型和程序是否正确实施或给出求解模型的误差、不确定性大小及使用范围,确认是要通过数值结果回答物理模型是否反映了真实客观世界或反映真实客观世界的可信程度.文章围绕爆轰流体力学模型,剖析了模型中不确定性因素,给出了影响模拟结果不确定性的关键因素清单,并对其开展了敏感度分析,确认了模型的适应性.     

抽运-检测型非线性磁光旋转铷原子磁力仪的研究

报道了一种抽运-检测型的非线性磁光旋转铷原子磁力仪.其原理是线偏振光通过处于外磁场环境中被极化的原子介质后,由于原子对线偏振光中左、右圆偏成分不同的吸收和色散,导致光的偏振方向会产生与磁场相关的转动.分析了该磁力仪的工作原理,并测试了它对不同磁场大小的响应.测试结果表明,磁力仪测量范围为100—100000 nT,极限灵敏度为0.2 p T/Hz~(1/2),磁场分辨率为0.1 p T.进一步研究了不同磁场下原子系综极化态的横向弛豫时间,讨论了原子磁力仪高磁场采样率的获得方法.本文的原子磁力仪在5000—100000 n T的磁场测量范围内磁场采样率可实现1—1000 Hz范围内可调,能够测量低频的微弱交变磁场.本文的研究内容为大磁场测量范围、高灵敏度、高磁场采样率的原子磁力仪研制提供了重要参考.     

抽样法与灵敏度法k_(eff)不确定度量化

核反应堆的中子学模拟计算中,核数据的不确定度导致的积分量计算结果的不确定度,通常采用基于微扰理论的灵敏度与不确定度分析方法 (简称灵敏度法)量化.灵敏度分析法原则上只适用于线性模型,且一般输运计算程序难以直接进行灵敏度分析.而抽样法直接抽样核数据输入中子学计算程序进行计算,通过对计算结果的统计分析评估计算量的不确定度.抽样法易于实现、计算精确、且适用性强.在灵敏度分析与不确定度量化程序SURE中,增加了抽样法不确定度的量化功能.为将抽样法不确定度量化应用于复杂问题的模拟计算,需对其进行细致的考核.为此,选取简单的临界基准实验模型,分别采用灵敏度分析法和抽样法进行不确定度量化,得到了各核素各反应道核数据导致的k_(eff)计算不确定度.对比显示,两种方法的不确定度计算结果有很好的符合,验证了SURE程序抽样法功能的正确性.抽样法计算的k_(eff)符合正态分布,说明在一般核数据的不确定度范围内,k_(eff)与核数据近似成线性关系,利用灵敏度分析法评估k_(eff)计算值的不确定度是适用的.